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We consider the effects of a perpendicular magnetic field on the d-density wave order and conclude 
that if the pseudogap phase in the cuprates is due to this order, then it is highly insensitive to the 
magnetic field in the underdoped regime, while its sensitivity increases as the gap vanishes in the 
overdoped regime. This appears to be consistent with the available experiments and can be tested 
further in neutron scattering experiments. We also investigate the nature of the de Haas- van Alphen 
effect in the ordered state and discuss the possibility of observing it. 



PACS numbers; 

Recently it was arguedEB that the observed pseudogap 
in the cuprate superconductors is due to the develop- 
ment of a second type of order, the cZ-density wave order 
(DDW), which is a particle- hole condensate with an in- 
ternal "angular momentum" 2. The notion of this broken 
symmetry has allowed us to understand an array of ex- 
periments, although a priori there is no sense in which 
such a state is close to, or adiabatically continuable to, 
a superconductor, or a Fermi liquid, which are the two 
most prominent states of matter. 

From considerations of the simplest Hartree-Fock 
Hamiltonian that captures the broken symmetry of the 
DDW and the electromagnetic gauge invariance, we 
study the dependence of the DDW order parameter as 
a function of the applied magnetic field perpendicular to 
CuO-planes. We find that the effect of the magnetic field 
on the pseudogap is rather weak, if the DDW order is 
well-developed. For example, to destroy a DDW gap of 
magnitude of order 20 meV would require fields of order 
1000 T, or larger. Moreover, this gap is highly insensitive 
to fields of order 25 T. This is, of course, a characteris- 
tic of the underdoped regime. In the overdoped regime, 
where the DDW gap drops rapidly to zero, the depen- 
dence on the magnetic field could be substantial. 

We suggest that the recent neutron scattering exper- 
iments, which report a tantalizing evidence of DDW 
orderB in the form of a resolution-limited elastic Bragg 
peak at the in-plane wave vector (7r/a,7r/a), where a 
is the lattice spacing, should be carefully examined for 
its dependence on a perpendicular magnetic field. The 
prediction is that the Bragg scattering intensity will be 
hardly affected, if the zero temperature DDW gap is 
substantialD 

The effect of magnetic field on a superconductor is well 
known. Since the Cooper pairs involve time reversed 
mates, the effect of any time reversal breaking perturba- 
tion is very important. In contrast, the particle- hole pairs 
of the same spin orientation form the DDW condensate, 
and, as a consequence, such a condensate is inherently 
less affected by a magnetic field. 

The contrasting response of a d-wave superconductor 
(DSC) and a DDW to applied magnetic field is further il- 
luminated by examining the nodal quasiparticles of these 
two systems. The nodal quasiparticles of a DSC do not 



form Landau levels,Ll while those of a DDW state do.El 
This implies that, in principle, it is possible to observe 
de Haas- van Alphen (dHvA) oscillations in a DDW state, 
but not in a d-wave superconductor. The reason is that 
quasiparticles in a superconductor do not couple simply 
to the vector potential A corresponding to the magnetic 
field B = VxA, but to the supercurrent {'V(p—2eA/hc), 
where ip is the phase of the superconducting order param- 
eter, e is the electronic charge, and c is the velocity of 
light. We calculate dHvA oscillations of the magnetiza- 
tion in the DDW state and find that although AB/B 
could be of reasonable magnitude, the amplitude of the 
magnetization oscillations is weak, generically of order 
O.Ol/is. Thus, observation of these oscillations may be 
difficult. 

Consider the mean field Hamiltonian of a d-density 
wave subjected to a magnetic field applied perpendic- 
ular to the plane. With the choice of the Landau gauge, 
the Peierls substituted Hamiltonian is written as: 

H = — ^ ^ ^ C^i^n^mn ~ t ^ ^ (c^ „Cm-|-l,Ti + h.c) 
m.n m,n 

^ <^(e2--"c+,„c„,,„+i+h. c.) 

7n.n 

- $](-l)"+"A (-^c+,„c„+l,„ + h.c) 

+ ^(-l)"+"A(-*e2— c+,„c™,„+i-|-h. c.) 

(1) 

Here m and n are the site labels of a square lattice in 
the X- and y-directions; a = 4>/4>o is the magnetic flux 
through each plaquette of the square lattice in terms of 
the flux quantum (/jq — hc/e, and Cm,n is the fermion 
destruction operator at the site {m,n). Here on, the 
hopping matrix element t will be set to unity. The DDW 
gap satisfies the self-consistency condition: 

A = -z(-l)™+"l/(c+,„c™+i,„ - h. c), (2) 

where the angular brackets denote the groundstate aver- 
age. Due to the conservation of current, ensured by the 
gauge invariance of the Hamiltonian d^), the DDW gap 
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could as well be: 

A = *(-l)"+"y(e2-"'"c+,„c,„,„+i - h. c). (3) 

Since a uniform applied magnetic field cannot break 
translational invariance, all physical quantities must be 
translationally invariant. Here V is an energy parameter 
that controls the strength of DDW pairing. 

We rewrite the Hamiltonian in terms of the Fourier 
transformed variables, given by (with lattice constant set 
to unity) 

c™,„ = ^ E (4) 
kk 

where N is the total number of sites on the lattice, to get 

- A ^ («e'''''cL+.+2.a,fc„+.Cfc^,fe« +h. c.) . (5) 



and the matrix i7k is 



One sees immediately that the states {kx,ky), [k^ + 
27ra, ky), {k^ + 4™, ky), . . . , (fc^; + tt, fc^ + tt), {kx + iT + 
2na, ky + tt), (/cj; +7r + 47ra, fcj, + 7r), ... are coupled. For 
rational flux, a — p/q, where p and q are relative primes, 
the set of coupled states is finite, containing 2q elements- 
At this point, the connection to the Hofstadter problemEl 
becomes particularly obvious. The Hamiltonian can be 
written as the following quadratic form: 



(6) 



where k^ G [0, 27ra) and ky e [0, tt). The row vector 
is given by 



— {^k^,ky I Cfe^-)-27rQ,fey i (^k^+iTTa,ky j • ■ ■ i (^k!c+2TTa{q—l),ky 7 
Cfc^ +7r,A,'y +7r ; ^kx-\-'7T-\-27Ta.ky+7T : ^kx+TT-i-ATTa.ky+TT 1 • • • J 
^kx -|-7r-|-27rQ(q — 1) jfcy -{-it ) : 



-2cosfc2; — e 

■^ky -2cos(fc2; + 27ra) 



-e 




2zAcosA:a; —iAe 
-iAe'-'^y 2iA cos(fc^ + 27ra) 
-iAe^^y 





-2 cos(/ca; + Ana) 


— iAe"''^" 



_^iky 



-2iA cos fc. 




zos{kx -\ 



zAe*'^" 



iAe*'=« -2iAcos(fc:^, + 27ra) iAe"*'=« 



2 cos /C;l 



2 cos(fca; + 27ra) e 







(7) 



The gap satisfies the equation 



obtained as 



<+^,fc .^Cfc^,fcJ +h. c. 



V 

N 



(8) 



For each value of {kx,ky), the Hamiltonian was diag- 
onalized, giving rise to a new set of fermion opera- 



,2q. Then, A 



tors {d^}, i — I, 

■ ■ ■ 1 ^k^-\-2TTa(q—l).kyi Afc^4.7r,A;y+7r j 



Afe 



-\'7T-\-27ra^ky +7r ; 



Afc^+^+2^a(5_i),fc^+^ were computed and so was the en- 
ergy spectrum. Several values of k^ and ky were chosen 
in the sub-Brillouin zone (sBZ) and finally the gap was 



2qNxNy 



(9) 



kesBZ 



where N — {qNa:){2Ny). The whole process was repeated 
until A converged. 

Typical magnetic fields in laboratories are of order 
20 T, equivalent to (l/1225)0o for the cuprates with lat- 
tice spacing of 4 A. In our computation, we were able 
to diagonalize matrices of size 4000 x 4000, or g = 2000 
corresponding to the smallest field of 12.5 T. Since the 
gap is so insensitive to the applied field, going to smaller 
fields appeared unnecessary. 

In Fig. we show the dependence of the DDW gap 
versus magnetic field at absolute zero for ^ = — O.lt as a 
function of V . We observe a striking insensitivity of the 
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FIG. 1: The DDW gap magnitude A as a function of the 
applied flux (j}/4>0 for various values of the coupling strength 
V/t (O : 0.5; □ : 0.75; O : 1.0; A : 1.25; < : 1.5; V : 1-75; 
> : 2.0) for = -Q.lt. 

gap with respect to the magnetic field. The gap only col- 
lapses, through a first-order transition at a flux strength 
4>cl 'Poi corresponding to a large laboratory magnetic field 
((^0 ~ 25, OOOT for cuprates). The gaps at zero field was 
calculated separately, and the results for finite field were 
seen to converge to them as ^ 0. It is important to 
note that because the chemical potential is finite, the 
DDW transition does not take place at F = O"*" , but at a 
finite threshold value Vc. This is clearly seen in Fig. ||. In 
the same figure, we see that (pc/cf'o drops quite rapidly to 
zero as Vc is approached. Thus, although ipc corresponds 
to quite large magnetic fields when the DDW gap is large 
in the underdoped regime, it could be quite small when 
the gap is small in the overdoped regime. For a gap of 
order 20 meV, the critical magnetic field is of order 1250 
T. 

In order to test the correctness of our computations, 
we also calculated the groundstate energy density at half- 
filling (/i = 0). We observe an increase in the groundstate 
energy as the magnetic field is applied; this increase in- 
dicates a diamagnetism in the DDW svstem. This result 
has, in fact, been analytically derivedQ using the nodal 
fermi gas formalism of the DDW state for fi ~ 0. In par- 
ticular, it was shown that the DDW groundstate energy, 
expressed in terms of lattice parameters, behaves as 

E,M)^Eas«»^,(^A,M(^f, ,10, 

where C(a;) is Riemann's zeta function, A and W be- 



ing the system area and a cutoff energy. The above for- 
mula is expected to hold only for small fields such that 
hujc — heB/m*c << W. This appears to be a surpris- 
ingly robust result, as shown in Fig. ^ especially because 
we made no approximations involving nodal Fermions. 
The fit at the smallest field gave a sensible W = 1.6t. 
The quantity W was interpreted in Ref. ^ as the cutoff 
energy below which the Dirac fermion description holds 
for the DDW state. 



10-' 



O 

o 

LU 

o 
m 

" 10-^ 



10"^ 10"^ 10"' 

FIG. 3: The groundstate energy as a function of the applied 
magnetic field at jj. — 0. 

As a consequence of the nonanalytic field dependence 
of the energy at /it = 0, the diamagnetic susceptibility 
diverges at zero field: x B'^/"^. Obviously, this di- 
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vergence poses a serious question about the stability of 
the DDW state at /x = 0. It is clear, however, that this 
divergence will be cut off for a finite chemical potential, 
fi ^ 0. Similarly, the divergence will also be removed 
at finite temperatures even for fi = 0. Moreover, three 
dimensional coupling between iJae layers will also modify 
the result for the susceptibility.B It is important to stress 
that the divergence in the susceptibility is intimately re- 
lated to the appearance of zero-energy modes, which is a 
direct consequence of the nodal points and the doubling 
of the Dirac fermions. 

It has recently been proposed in Ref. || that a finite 
magnetic field may induce an additional dxy gap. Such 
a gap can efficiently remove the nodes and hence the di- 
vergence in the diamagnetic susceptibility. However, the 
small dxy gap will not be manifest at a finite temper- 
ature and the dominant cutoff will be provided by the 
temperature or the chemical potential. Therefore our 
computation for the DDW gap retains its relevance to 
experimental situations. 

As an interesting contrast to DSC, an essential feature 
of the Dirac fermion picture of the EtDW is th e formation 
of Landau levels in a magnetic field u £"„ = y/n{HuJc)i!^o, 
where Aq is the amplitude of the gap, and ujc is the cy- 
clotron frequency, eB/'m*c, m* being the effective mass 
corresponding to the nodal region. As a result, the DDW 
system away from half-filling is expected to exhibit de 
Haas-van Alfven (dHvA) oscillations in the magnetiza- 
tion. 

For a free electron system at chemical potential /z, 
where the energy spacing is huc, the oscillation frequency 
satisfies 



■0.04 



(11) 



Taking B = 1 T, ^ « 0.25 eV, one would get AB « 
5 G. For nodal fermi gas at chemical potential fj,, the 
oscillation frequency is 

To make an estimate, we can choose B — IT. It is more 
problematic to estimate /i in the underdoped regime in 
which the DDW state is likely to be found, and in which 
the typical non-interacting band picture may not hold. 
There is some CYidence, however, that /i w — te^, where 
X is the doping.B Choosing x ~ 0.15, Aq of the order 
of the pseudogap ^ 20 meV, and t « 0.25 eV, results 
in AB « 0.05 T, which is clearly within the detectable 
range of dHvA oscillation experiments. 

The above estimate for AB is not sufficient for the de- 
tectability of dHvA effect. We must also estimate the 
magnitude of the magnetization oscillations. For this 
purpose, it is convenient to consider energy spectrum of 
the nodal fermions of the DDW state. The oscillation 
in energy, and thus of the magnetization, occurs when- 
ever a Landau level passes through the Fermi level. For 
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FIG. 4: The de Haas-van Alphen oscillations of the magneti- 
zation in the DDW state. 



hole-doped cuprates, with a finite chemical potential /z, 
the formula (^) for the groundstate energy must be aug- 
mented by the following additional term: 



eB 
he 



(13) 



T!=0 



where [rifj] is the largest integer smaller than ij?/ (SwcAq). 
In other words, is the number of filled Landau lev- 
els. With the above choice of parameters, the magneti- 
zation per unit cell, in units of Bohr magneton, /zs, is 
shown in Fig. ^. In plotting this figure, we have removed 
the spurious 5-function singularities arising from taking 
derivatives of the step functions. At finite temperature, 
the sharpness of the oscillations will be smoothed out, 
but the oscillation amplitude will still roughly be of or- 
der « 0.01 /is, which is quite weak. As is evident from 
Fig. |[ as the magnetic field swept between 0.5 T and 1 
T, the system undergoes about 14 oscillations. 

There are already some experiments that appear to 
shed light on the effect of magnetic field on the pseu- 
dogap. Nuclear magnetic resonance experiments in un- 
derdoped YBa2Cu408 clearly support the view that 
the pseudogap is hardly j-aJfected by an applied per- 
pendicular magnetic fieM,E3 and so does experimets in 
nearly optimally dopedllil YBa2Cu307_5. In contrast, 
nuclear magnetic resonance experiments in overdoped 
TlSr2CaCu206.8 provides contrasting evidence that in 
this pigime the pseudogap is sensitive to the magnetic 
field.ll3 At first sight, the experiments in YBa2Cu307_5 
appear to be coutipj/ersial because there are also reports 
to the contrary,t3o but as was pointed out by the au- 
thors of Ref. [ij, the later reports by the authors of 
Ref. |l3|,^ seem to be in agreement with Ref. |ll|. There- 
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fore, our theoretical results would be consistent with ex- 
periments, if the pseudogap is indeed caused by DDW. 
We hope that our work will catalyse further experiments 
on this important question. On the theoretical front, it 
will be necessary to go beyond the Hartrce-Fock approx- 
imation, especially at finite temperatures. To go beyond 
Hartree-Fock approximations, we must recognize the fact 
that the DDW gap, which characterizes the orbital cur- 



rents need not be in the perfect staggered pattern, but 
can be flipped as long-as the current conservations at the 
vertices are satisfied.Ej We hope to report on this issue 
in the future. 
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